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Quantities in columns 2, 3, 4 and 5 are from de Vaucouleurs et al. (1991) (RC3) 

References for Table 1. 
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mann (1981); (4) Peterson (1978); (5) Bottinelli et al. (1984); (6) Kennicutt 
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Table 3. Star formation Histories of GEHRs 
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ABSTRACT 

A study of star formation is carried out on 35 giant extragalactic Hll regions 
(GEHRs) in seven galaxies using optical photometric data in BVR broad bands and 
in the emission line of Ha. Interstellar extinction, metallicity and nebular contributions 
to the broad bands are estimated using spectroscopic data on these objects. Dimension- 
less diagrams involving B — V and V — R colors and the flux ratio of Balmer line to B 
band continuum are used to study star formation. The cluster colors indicate reduced ex- 
tinction towards stellar continuum compared to the values derived from Balmer lines for 
the ionized gas. The frequency of detection of classical young {t < 3 Myr) regions with 
only one burst of star formation is found to be low as compared to young regions with 
an accompanying population rich in red supergiants from a previous burst {t ~10 Myr). 
Reduced extinction towards cluster stars, destruction of ionizing photons and the exis- 
tence of older population, often spatially unresolvable from the younger population, all 
conspire to make the observed Balmer line equivalent widths low in a majority of the 
GEHRs. A scenario of star formation is suggested which explains many of the observed 
properties of GEHRs, including the core-halo structure, reduced extinction for the ra- 
diation from stars as compared to that from the nebular gas, non-detection of young 
single burst regions and the co-existence of two populations of different ages. 
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1. Introduction 

Giant extragalactic H ii regions (GEHRs) are sites of 
current star formation activity in galactic disks. A 
detailed investigation of star formation in these rel- 
atively simpler systems is important towards under- 
standing star formation in more complex star-forming 
systems such as blue compact galaxies (BCDs or Hii 
galaxies), and starburst systems. GEHRs in an indi- 
vidual galaxy are known to differ from one another in 
many of the observational properties, but still have 
several properties in common. The common proper- 
ties such as the constancy of luminosity and diameter 
of the brightest (also largest) Hii regions have been 
used earlier for estimating the distances to galaxies 
(Sandage & Tammann 1974, Kennicutt 1979). It is 
important to investigate the reasons for similarities 
as well as differences between different H ii regions in 
order to make use of GEHRs as reliable distance in- 
dicators or tracers of star formation rates in galaxies. 

Initial Mass Function (IMF), age, history, metal- 
licity, extinction, mass or luminosity are among the 
most important parameters determining the proper- 
ties of GEHRs. In recent years, there have been sev- 
eral attempts to investigate the IMF parameters un- 
der a variety of conditions (see Scalo 1986). Yet, the 
IMF slope of a = 2.35 as obtained by Salpeter (1955) 
remains the most accepted value. Humphreys & 
McElory (1984) have favoured a much steeper IMF for 
high mass stars in the Galaxy, while Parker & Gar- 
many (1993) obtained an IMF slope close to Salpeter's 
value for the stars in 30 Doradus, both studies being 
based on star counts. Investigations based on the in- 
tegrated properties of GEHRs have so far remained 
inconclusive on the exact value of IMF slope and the 
lower and upper cutoff masses. There are some in- 
dications that IMF may be flatter or deficient in low 
mass stars in active starburst regions (Scalo 1987). 

In this study we aim to investigate the star forma- 
tion histories of GEHRs, taking into accoimt the pos- 
sible IMF differences from one region to another. It is 
however important to take into account the effects of 
(i) interstellar extinction, (ii) metallicity and (iii) neb- 
ular contributions on the observable quantities, be- 
fore these quantities can be compared with the syn- 
thetic models. Interstellar extinction is convention- 
ally derived using the observed ratios of Balmer lines 
(see e.g. Shields 1990). Balmer decrements as well as 
mctallicities and nebular contributions can be derived 
from spectroscopic data on individual regions. With 



the above problems in mind, we use the BVKHa 
data of only those GEHRs from Mayya (1994) for 
which additional spectroscopy is available in the lit- 
erature. The data, after applying corrections for ex- 
tinction and nebular contribution, are compared with 
the evolutionary population synthesis model to inves- 
tigate the properties of star formation in the sample 
regions. Analyses were performed on dimensionless 
diagrams, thus eliminating the possible errors due to 
the uncertainties in the estimation of distance to the 
source and/or the mass of the star-forming region. We 
find that the extinction experienced by the embedded 
cluster stars is lower compared to the values derived 
from the Balmer decrement. These results are of spe- 
cial interest, especially after a similar conclusion has 
been reached by the independent study of Calzetti et 
al. (1994) using ultraviolet (UV) and optical spectra 
of central regions of galaxies. Our study also reveals a 
slightly complex history of star formation in GEHRs. 

In sec. 2, we describe the sample and the method 
followed in extracting the pure stellar continuum from 
the observed BVR data. In sec. 3. the properties of 
the embedded cluster are compared with population 
synthesis models. Our results are discussed in the 
context of recent results in sec. 4. Physical mod- 
els for star formation in GEHRs and the distribution 
of cluster stars, gas and dust are proposed in sec. 5. 
Sec. 6 summarizes the main results of this study. 

2. Sources of Data and Preliminary Treat- 
ment 

The GEHR data for this study are taken from the 

BVRHa photometric catalogue of Mayya (1994, Pa- 
per I henceforth). The catalogue contains CCD aper- 
ture photometric data on nearly 180 regions in nine 
nearby galaxies. In these studies, special care was 
taken in ensuring that the contamination by the back- 
ground is minimum and hence the derived magnitudes 
are genuinely from the embedded cluster and the sur- 
rounding nebulosity. The usage of apertures rather 
than slits, as well as measurement of BVR continuum 
and Ha emission line fluxes, centered on the same 
pixel, makes these measurements ideally suited for 
the investigation of star formation properties. Spec- 
troscopic data are available for about one fourth of 
these regions (~45). The sensitiveness of the back- 
ground values on the catalogued colors is indicated 
by a photometric quality index in the catalogue. Re- 
gions having photometric quality index 4 (background 
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uncertainties leading to errors in exc;css of 0.3 mag in 
B — V and V — R colors: see Appendix in Paper I for 
a discussion on background subtraction errors in H II 
region photometry) arc dropped to eliminate back- 
ground dominated regions from the sample. The es- 
timated rms errors for the entire sample in the cata- 
logue arc < 0.15 mag in colors and < 0.05 dex in the 
ratio of H/3 to blue band luminosity {^)- The re- 
gions selected for the study here are relatively brighter 
and hence the errors are expected to be smaller. We 
estimate the errors to be a factor of two less as com- 
pared to that for the entire sample. Finally, we are 
left with 5 regions in NGC 1365, 3 in NGC 2903, 14 in 
NGC 2997, 2 in NGC 3351, 3 in NGC 4303 in addition 
to the giant Hii complex NGC 2363 in NGC 2366 and 
the starburst nucleus of NGC 5253. Measurements 
using bigger apertures, which include more than one 
emitting knot, are also included around six regions in 
NGC 2997. 

The spectroscopic data involving [O iii] and Balmer 
lines are used to estimate the abundance of oxy- 
gen, electron temperature (Tg) and visual extinction 
(^^(Bal)) using a semi-empirical approach (McCall 
et al. 1985). All the observed quantities have to be 
corrected for interstellar extinction before comparing 
with the models. Total extinction towards extragalac- 
tic Hll regions consists of three parts — (i) extinc- 
tion from the interstellar medium (ISM) of our galaxy 
(At, (gal)), (ii) extinction from the ISM of the parent 
galaxy outside the Hli region (^,,(ext)) and (iii) ex- 
tinction due to the dust mixed with the gas within the 
Hll region (yl„(int)). A„(gal) values for the program 
galaxies are given in Table 1 for each of the program 
galaxies, as tabulated by de Vaucouleurs et al. (1991). 
Other quantities in the table are distances to the par- 
ent galaxies, number of regions used in the present 
study and the sources of spectroscopic data and dis- 
tances. We denote the visual extinction as derived 
from the Balmer decrement by Ay, in the remaining 
part of the text, though it contains contribution from 
all the three components. We start our study using 
Ay to deredden the cluster and nebular quantities, as- 
suming that the spcctroscopically derived values are 
valid for both the cluster continuum and the nebular 
radiation integrated over the apertures. In the later 
part of this section, we suggest alternative prescrip- 
tions for extinction corrections that are necessary to 
explain the present data set. 

The usage of Ay to deredden quantities related to 
embedded clusters is appropriate only if the amount 



of dust residing internal to H II regions is negligi- 
ble. The large spread in extinction values for regions 
within a single galaxy implies non-uniform distribu- 
tion of dust within the parent galaxy, either within 
or outside the Hii region. There have been stud- 
ies investigating the relative contribution of the dust 
mixed with the ionized gas, known as 'internal extinc- 
tion', to the total extinction from the parent galaxy 
(Lequeux et al. 1981). These studies show the inter- 
nal extinction to be non-negligible at least for some 
regions. The extinction values derived using the radio 
thermal continuum is statistically known to be higher 
than that derived using the Balmer decrement (Israel 
& Kennicutt 1980). Caplan & Deharveng (1986) have 
studied the effect of different dust configurations, in- 
cluding the dust internal to the Hii regions to ex- 
plain the observed discrepancy between the radio and 
optically derived extinction values. Internal redden- 
ing is found to be significant for some regions in the 
Large Magellanic Cloud (LMC). We find the extinc- 
tion properties of embedded cluster stars to be dif- 
ferent from that of the ionized gas and hence suggest 
an alternate prescription for extinction correction in 
GEHRs. 

The Ha band fluxes have contribution also from 
[Nil] lines at A 6548, 6583 A. The observed Ha-h[Nii] 
fluxes are first corrected for interstellar extinction 
using Ay and the galactic extinction curve (Seaton 
1979). Ha fluxes are computed from the observed 
Ha-|-[Nii] fluxes using the spcctroscopically deter- 
mined [Nil] /Ha ratios. Lyman continuum luminosi- 
ties are derived from extinction corrected Ha lumi- 
nosities using 



N. 



^ j = 7.32 X 10" 



Ha 



(1) 



phs " erg s 

which assumes an ionization bounded case B nebula 
(Osterbrock 1989). The expected H/3 luminosity {(j)) 
is computed from this using 



4.78 X 10 



-13 



Nr 



(2) 



ergs " phs \ 

The treatment of broad-band magnitudes is slightly 
different. First, an estimation of the contribution 
from the nebula (line -|- continuum) within each band 
is made based on the derived Tg and extinction cor- 
rected Ha luminosity. The strengths of all bright 
lines relative to H/3 line, are taken from McCall et al. 
(1985). Free-free, free-bound and 2-photon emission 
mechanisms are considered in computing the nebular 
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continuum. Details of these computations are pre- 
sented in Mayya (1995). For a majority of the regions, 
the nebular corrections amount to < 5% in B and V 
bands and 15 25% in the R band, with emission lines 
being the main contributors. The estimated nebular 
contributions to B and V bands exceed that in the 
R band in the two high excitation regions (NGC 2363 
and 5253), due to the presence of strong [Oili] lines. 
However, these estimates heavily depend on the as- 
sumed shape of the passband, because of the location 
of the [O III] lines in the steeply falling part of the B 
and V band responses. We find that even a 50A shift 
of the V band response curve to longer wavelengths 
changes the B — V and V — R colors by as much as 
0.3 mag in NGC 2363. This amount of shift of the ef- 
fective wavelength is not unusual and hence gaseous 
corrected cluster colors using broad bands are bound 
to have large errors for high excitation regions. Use of 
narrower bands situated at emission line free regions 
would have eliminated this source of uncertainty. 

The estimated nebular luminosities are subtracted 
from extinction corrected broad band luminosities to 
obtain pure cluster quantities. These quantities are 
reddened back to obtain the observationally expected 
values without any contamination from the surround- 
ing nebula. We denote these quantities by the sub- 
script *, while the observed quantities are denoted by 
the subscript 'obs' {e.g. {B - V)ohs, (B - V)^ etc.). 
The resultant quantities are corrected for the galac- 
tic and total extinction using Ay (gal) and Ay and are 
denoted by suffixes 1 and 2 respectively. The Galac- 
tic extinction curve as tabulated by Seaton (1979) is 
used in these computations. The process is summa- 
rized below. 

(B-y)obs gassubt ^ (B-V)*, (3) 

{B -V). A,(gal) ^ (B - V)i , (4) 

{B -V% Av(Bal) ^ (B - V)2 , (5) 

Table 2 contains the results for individual H ii regions. 
Column 1 contains identification numbers for GEHRs 
following Paper I. Oxygen abundances in column 2 
are derived from ratio using the semi- 

empirical calibration of McCall et al. (1985). Note 
that the sample regions have metallicities in the range 
2-4^0, with the exception of the two high excitation 
regions mentioned above. McCall et al. also provide 
relation between the oxygen abundance and T^, , which 
is used to derive given in column 3. Column 4 lists 
the visual extinction based on the Balmer decrement. 



Oxygen abundance and visual extinction are based 
on the sources of spectroscopic data listed in Tab. 1. 
While deriving the visual extinction, the underlying 
absorption in Balmer lines, which at the most con- 
tributes about 2 A in H/3 equivalent width (McCall 
et al. 1985), has been accounted for. The K magni- 
tude, and {B — V)^, and {V — i?), colors in columns 5, 
6 and 7 respectively, are cluster quantities after sub- 
tracting the gaseous contribution. Column 8 contains 
the log(-^)» = \og{(t>) — log(LB)^ values. Note that 
Ay (Bal) has been used in obtaining the extinction cor- 
rected (j) from observations. Lyman continuum lumi- 
nosity is given in column 13. Remaining columns in 
the table will be explained in the next section. 

3. Comparison with Models 

We plot in Fig. 1 the values of log(-j^)i against 
{B — V)i, which are corrected only for galactic extinc- 
tion. Fig. 2 shows a similar plot, but with {V — R)i 
replacing (B — V)i. Hii regions from different galax- 
ies are shown by different symbols, as indicated at the 
right hand side of the diagram. The small and large 
aperture measurements are distinguished in the case 
of NGC 2997. Rms errors on the plotted quantities 
are indicated by the cross at the bottom-left corner 
of the diagram. For the two high excitation regions, 
NGC 2363 and 5253, there is an additional source of 
error, arising out of the subtraction of nebular con- 
tribution, as discussed in the previous section. The 
arrow on the right side of the figure is the reddening 
vector, with its length corresponding to 1 magnitude 
of visual extinction. 

The observed quantities are compared with the 
evolutionary population synthesis model, which is 
briefly described below. The computational details 
and results for different star formation scenarios can 
be found in Mayya (1995: Paper II henceforth). The 
evolutionary population synthesis model is based on 
the input stellar evolutionary data of Schaller et al. 
(1992), and stellar atmosphere data of Kurucz (1992). 
The cluster is defined by an initial mass function, age, 
metallicity and the total mass in stars. The evo- 
lution of the cluster is performed at every 0.1 Myr 
interval, which is adequate to sample even the fast 
evolutionary phases of massive stars. Optical and 
near infrared colors are computed and tabulated in 
Paper II for three different scenarios of star forma- 
tion at solar metallicity. As noted earlier, GEHRs 
in the present sample have metallicities between 2- 
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4Zq, with the exception of NGC2363 and 5253 
and ^Zq respectively). Recently, Cervifio & Mas- 
Hesse (1994) and Leitherer & Heckman (1995), have 
shown the metallicity differences to affect the clus- 
ter evolution significantly. Hence we extended the 
model computations of Paper II to Z= O.IZq and 
Z= 2Zq based on the stellar evolutionary tracks of 
Schaller et al. (1992) and Schaerer et al. (1993) re- 
spectively. Metal-poor models are particularly impor- 
tant for the study the star formation in blue compact 
galaxies. Note that the metallicity of the most metal- 
poor galaxy, IZwl8 is around a factor of 4 smaller 
than our metal-poor model. The evolutionary tracks 
for an Instantaneous Burst (IB) of star formation at 
the three metallicities are shown in Figs. 1 and 2. 
The long-dashed curve represents the solar metallic- 
ity model, while the short-dashed and the solid curves 
represent O.IZq and 2Zq models respectively. The 
clusters arc evolved upto 6.5 Myr, with dots placed 
on the curves at every 0.5 Myr. The evolutionary 
phase corresponding to 4 Myr is indicated by the let- 
ter 4 on each of the curves. The IMF chosen is close 
to Salpeter's: m; = 1 Mq, niu = 60 Mq and a = 2.5. 

Two phases can be distinguished during the early 
(age<10 Myr) evolution of an IB. The first phase is 
controlled by the massive main sequence stars while 
the second by the supergiants, in particular the red 
supergiants. The transition from the first phase to 
second phase takes place around 4-5 Myr, with metal- 
licity playing a major role during both the phases. 
The main sequence life time of the stars decreases 
with the increasing metallicity resulting in a faster 
drop of in metal-rich regions during the first 
phase. Higher metallicity favours the formation of 
red supergiants (RSGs) and Wolf-Rayet stars in the 
cluster and hence a given burst of star formation pro- 
duces more of these stars in a metal-rich environment. 
Additionally the metal-rich RSGs are cooler than the 
metal-poor RSGs, resulting in redder colors during 
the red supergiant phase of metal-rich clusters. The 
increased optical luminosity due to RSGs in metal- 
rich regions also leads to a decrease of The 

metal-poor locus turns back at 5 Myr (-^ increases 
and colors become bluer), which can be imderstood in 
terms of the absence of long lived cooler supergiants in 
the cluster. Thus, the metal-rich clusters have lower 
values of throughout their evolution, and redder 
colors after around 4 Myr of evolution in comparison 
to metal-poor regions. A majority of the regions in 
our GEHR sample are metal rich and hence we use 



Z=2Zq models for comparison with our data. 

For m„ other than 60 (stellar masses are expressed 
in solar units henceforth) the evolutionary track is 
similar to the one shown, except that the starting (age 
= 0) value is different. For m„ > 60, it starts at a 
higher value than shown and the converse is true for 
lower masses. At a = 2.5, the B — V color reddens 
by 0.04 mag as the upper cut-off is decreased from 
120 to 3OM0, whereas log(i^) decreases by 0.70 dex 
for the same range of masses. As the more massive 
stars die, the evolution follows that of a cluster with 
a lower m„. 

It can be noticed that the evolutionary track forms 

an envelope around the observed points. All the ob- 
served points can be brought into the locus expected 
from the Z= 2Zq model, by shifting the points along 
the reddening vector by amoimts Ay =0-1.4 mag. 
The giant H 11 complex NGC 2363 has several proper- 
ties which distinguishes it from the rest of the sample, 
and hence we treat this region separately in the dis- 
cussions to follow. The distinguishing properties are 
its low metallicity, negligible extinction (Kennicutt et 
al. 1980), and a large correction for the nebular con- 
tribution in B,V and R bands. 

In Figs 3 and 4, the observed quantities are plot- 
ted after dereddening by Ay, as derived from Balmcr 
decrement (quantities with suffix 2). Other details 
in these figures are similar to Figs 1 and 2, except 
that only Z= 2Zq evolutionary track is shown and 
the track is truncated at 6 Myr instead of 6.5 Myr. 
It is striking to note that the reddening correction 
has moved the points away from the models in both 
the axes. Hence the corrected B — V and V — R 
colors are too blue and (7^) too low compared to 
the values expected for clusters younger than around 
4 Myr. One can infer from Figs 1 and 2, that there are 
many regions for which very little extinction correc- 
tion is required, whereas the extinction values com- 
puted using the Balmer decrement have a median 
value close to 1.5 mag, with only 4 regions having 
values less than 0.75 mag. These two figures indicate 
that by dereddening the cluster related quantities by 
extinction derived from the Balmer decrement, we are 
over-correcting for extinction. In order to investigate 
this in more detail, we plot in Fig. 5, the dereddened 
{B — V)2 colors against the pure extra-galactic visual 
extinction Ay — Ay (gal). A trend of blue {B — V)2 
colors having higher extinction is clearly seen in the 
figure. A similar trend is also seen in {V — R)2 color, 
but to a lesser degree (not plotted). In fact both 
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{B — V)2 and {V — R)2 colors are bluer than the model 
for regions having Ay > 1.5 mag. The figure gives an 
important clue on the nature of interstellar extinction 
towards stellar continuum and ionized gas emission in 
GEHRs. The inescapable conclusion from Figs 3, 4 
and 5 is that the cluster continuum on an average ex- 
periences lesser extinction than the emission from the 
ionized gas. Such a situation is possible if most of 
the extragalactic extinction is due to the dust inter- 
nal to GEHRs and the radiation from the ionized gas 
is selectively absorbed by the dust. The ionized gas 
in GEHRs is known to be confined to clumps, with a 
volume filling factor as small as 0.01 (Shields 1990). 
The dust is also likely to be present in the same vol- 
ume. Radiation from the stars, which contribute to 
the optical and UV continuum, on the other hand 
might occupy less dustier regions, and hence expe- 
rience lesser amount of extinction. The segregation 
might be arising as a result of relative motion of the 
stars with respect to the gaseous material from which 
they have formed (Leisawitz & Hauser 1988). Based 
on this simple clumpy model, we formulate a different 
prescription to estimate the extinction of the stellar 
continuum which should be more accurate at least 
statistically. 

3.1. Estimation of Extinction towards Clus- 
ter Continuum 

The observed range of colors is too large to be ex- 
plained in terms of IMF variations. For an assumed 
slope of the IMF, the B -V color reddens by < 0.04 
mag as m„ decreases from 120 to 3OA^0. Color differ- 
ences are around 0.20 mag even for two extreme IMFs, 
one rich in massive stars (m„ = 120, a = 1.0), and 
the other rich in low mass stars (mi, = 30, a = 3.5). 
Thus, as a first approximation, we can neglect the 
dependence of colors on the IMF for clusters younger 
than 3.5 Myr. We assume B-V = -0.25 and V-R = 
—0.10 as typical values for the young clusters in our 
sample, corresponding to a = 2.5. and Z = 2Zq. The 
remaining part of the discussion pertains to the IMF 
parameters mi = 1,to„ = 60, a = 2.5, unless other- 
wise mentioned. The regions having high value 
in addition to blue B — V and V — R colors must be 
the youngest regions in the sample. We hence make 
independent estimate of visual extinction (A^c) to- 
wards stellar continuum on the basis that dereddened 
B — V and V — R colors cannot be bluer than —0.25 
and —0.10 respectively. The method is summarized 
by the following equations, which are based on the 



formulation of Calzetti et al. (1994). 

Consider an embedded cluster with intrinsic inten- 
sity 7° at a given wavelength. The intensity Ix leaving 
the parent galaxy is then given by. 



(6) 



where tx is the line of sight optical depth suffered by 
the radiation at wavelength A in the host galaxy. If 
Tfj and Tq, are the optical depths experienced by the 
Balmcr line photons H/3 and Ha respectively, then 
the differential optical depth between H/3 and Ha is 
given by. 



TfJ 



= ln 



f {Ha/H(i), 
\ 2.86 



(7) 



where (Ha/H/3)j is the ratio of intensities of Ha and 
H/3 lines after correcting for galactic extinction. Vi- 
sual extinction A^ tabulated in Tab. 2 is related to 



by, 



A,(Bal) - Av(gal) = — — , (8) 



or 



X/3 — Xq, 

= 0.354(A^(Bal) - Av(gal)) , (9) 



where Xp and Xa are tabulated by Scaton (1979) for 
the galactic extinction curve. Tg as defined above is 
the difference in extinction between H/3 and Ha emis- 
sion lines. Analogously we can define the extinction 
suffered by the continuum photons between H/3 and 
Ha wavelengths. We denote this quantity as de- 
fined as, 



Xj3 — Xa 

1.086 



E{B - V)i 



(10) 



where 



E{B - V)i = {B- VUs -{B- y)mod - A,(gal)/3.1. 

(11) 

We define. 



5= -r 



Ay {gal) 



Ay- Ay (gal) 



(12) 



Equations similar to the above are also formu- 
lated for the determination of g from the observed 
V — R colors. There are 11 regions in the sample for 
which both {B — V)2 and (V — R)2 are bluer than 
the bluest model colors. We derive a mean value of 
g = 0.70 ± 0.20, from both these colors. Physically g 
is controlled by the spatial distribution of gas, dust 
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and stars with a value of unity if all the three compo- 
nents are well mixed inside the volume of GEHR or 
the entire dust is situated far away from the emitting 
region. However there are several issues which come 
in the way of an accurate determination of g, some of 
which are intrinsic to the region of interest and hence 
unavoidable, while others are due to observational er- 
rors on colors. In estimating g, we assumed all re- 
gions whose dereddened colors become bluer than the 
bluest model colors to be young. However a moder- 
ately evolved region with a lower value of g than esti- 
mated is indistinguishable from a young region with g 
value closer to unity. Thus the actual value of g could 
be smaller than the estimated mean. In estimating g, 
we also assumed that the entire extra-galactic extinc- 
tion comes from within the star-forming complex. In 
reality, there may be some contribution from the ISM 
of the parent galaxy. Correction to this has the effect 
of decreasing the value of g. With all these consid- 
erations, we assume g = 0.50 to be typical for the 
whole sample of GEHRs in the subsequent analysis. 
This value is supported by a recent analysis involving 
UV continuum and Balmer lines of central regions 
of galaxies (Calzetti et al. 1994). Using this aver- 
age value, A^c is recomputed for all the regions using 
eq. 12. The resultant values are tabulated in column 
9 of Tab. 2. The cluster quantities are dereddened 
using these values of A^c and the resulting quanti- 
ties are identified by suffix 3. Nebular quantities and 
hence <f) are corrected for extinction using A^. 

(B-V), Avc^(B-V)3 (13) 

The resultant values of (-^)3, {B - V)3 and (V - 

R)3 are tabulated in Tab. 2. log(-^)3 is plotted 
against (B — V)3 and {V — R)3 in Figs 6 and 7 re- 
spectively. As in the earlier figures the IB model for 
the chosen IMF is shown with solid curve, with dots 
placed every 0.5 Myr. The reddest point on the curve 
corresponds to 6 Myr of evolution. With the new pre- 
scription for extinction correction, not only the colors, 
but also the log(i^) values are brought closer to the 
curve represented by the IB model. While a set of 
regions follow this track in both the diagrams within 
the observational errors, there are some regions which 
lie below and some others which lie to the right (red- 
ward) of the solid line. These regions require alterna- 
tive models, which is the main theme of the discussion 
in the remaining part of this section. 



3.2. Destruction and Escape of Ionizing Pho- 
tons 

There are statistically significant number of regions 
which have blue colors, indicative of young age, but 
have values lower than that expected for a clus- 
ter younger than 3 Myr. These regions seem to agree 
well with the model, if the log(-j^) is shifted by about 
—0.2 dex, as indicated by the dashed curve in Figs 6 
and 7. The possible sources, which can give rise to 
this shift are investigated below. We express the dis- 
crepancy in log(-j^) between observations and the 
model evolutionary track by a parameter / defined 
as, 

log(^)obs = log(^)mod(0 + log /, (14) 

where (;^)obs ={^)3 and log(;^)inod(i) is the IB 
model value at age t. For the IB model, changes in 
are mainly controlled by the changes in (j) and hence, 
physically, / is the ratio of the detectable number 
of Lyman continuum photons to that expected from 
the model and has a max;imum value of 1 when every 
Lyman continuum photon emitted by the cluster stars 
at the estimated age can be traced by the Balmer 
emission lines. We estimate ages between 2.5-4 Myr 
for the regions which require —0.2 dex shift. This 
value of shift implies a value of / = 0.65 or 35% of 
the ionizing photons could not be traced by the Ha 
emission. 

The following processes may be responsible for the 
missing ionizing photons: (i) the destruction of ion- 
izing photons by dust within the ionized complexes 
('?dusi), and (ii) the escape of ionizing photons from 
the nebula or an under-estimation of Ha flux, which 
traces the ionizing photon rate (?7esc)- The under- 
estimation may be caused by the total absorption of 
Ha photons by dust clumps or because a significant 
fraction of the total Ha luminosity is in low surface 
brightness halo component, which lies below the de- 
tection limit in most cases. Jjdust and ?7esc are related 
to / by the eq., 

/ = (1 - ??dust)(l - '?esc)- (15) 

Good estimates of %ust) Vesc are lacking both obser- 
vationally and theoretically. Smith et al. (1978) and 
Mezger (1978) have estimated ?7dust to be between 
0.7 and 0.2 for giant galactic Hll regions. Belfort 
et al. (1987) have assumed a useful mean value of 
jydust = 0.30, a value also used by many others (see 
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e.g. Mas-Hcssc & Kunth 1991). The observed value 
of / = 0.65, then corresponds to neghgible rjesc- The 
value of /, however, depends on the assumed IMF 
parameters with the dependence weakening beyond 
3 Myr. 

From these analyses, it may be concluded that the 
absorption of ionizing photons by dust is the main 
source of the discrepancy in (or Balmer line equiv- 
alent widths) observed in some regions, for the as- 
sumed IMF parameters. This is also supported by our 
conclusion based on cluster colors, that dust resides 
inside GEHRs. On the other hand a non-zero r^esc 
cannot be completely ruled out because of the uncer- 
tainties in m„, age and absorption fraction of ionizing 
photons by dust. The detection of extra-Hii region 
ionized gas in several irregular galaxies by Hunter & 
Gallagher (1992) and Hunter et al. (1993) has raised 
questions regarding the source of ionization of these 
diffuse regions. The Lyman continuum photons es- 
caping the Hii complexes may explain the observed 
ionization. 

3.3. Effect of Multiple Star Bursts 

The evolution of an instantaneous starburst with 
around 35% loss of the Lyman continuum photons, 

along with reduced effective extinction towards ra- 
diation from stars compared to that from the ionized 
gas, is found to explain the blueward envelope (within 
the limits of observational scatter) of regions in the 
log(i^) vs {B — V) and {V — R) diagram. However, 

there is considerable spread of colors at any given 
which cannot be explained by an instantaneous burst 
model. 

The departure of the colors from the IB model 
is in the redder direction, which can be understood 
in terms of a population rich in red supergiants, co- 
existing with the younger population. The two popu- 
lations of different ages are seen spatially resolved in 
30 Doradus (McGregor & Hyland 1981). The younger 
population provides the ionizing radiation whereas 
the red supergiant population contributes bulk of the 
optical continua. Stars in the mass range 40 15M0 
become red supergiants at ages 4-14 Myr a,t Z = 2Zq 
(Schaerer et al. 1993). Simultaneous existence of ion- 
izing stars as well as red supergiants implies extended 
duration of star formation either proceeding continu- 
ously or in the form of bursts of short durations. The 
range of computed quantities for continuous star for- 
mation is small compared to the range of observed 



quantities (Paper II), and hence the continuous star 
formation cannot explain present observations. On 
the other hand, 2-burst models of Paper II reproduce 
much wider range of observable values. We discuss 
these models in order to explain our observations. 

The composite model of two starbursts of equal 
strength with the younger one at 0, 3 and 4.5 Myr 
is shown in Figs 6 and 7 as three sequences of short- 
dashed lines. Both the bursts assume IMF slope of 
2.5 and m„ = 60 and Z = 2Zq. In each of these 
sequences, the age of the older burst changes from 
6-14 Myr, with beginnings of the sequences marked 
by the short vertical line followed by successive solid 
dots spaced every 0.5 Myr. The color for IB model 
becomes reddest at around 7 Myr corresponding to 
the RSG phase of stars around 25M0 mass, turning 
blue again at later stages of evolution. This is the 
cause of the blue-ward turnover of the 2-burst mod- 
els. The general trend of many observed regions lying 
redward of the IB model is well reproduced by the 2- 
burst models, with younger burst ages between and 
4.5 Myr and the older burst in the supergiant-rich 
phase (age ^ 10 Myr). 

For the regions with two bursts of star formation, 
the observed values of log(-^) can be reproduced by 
models without the need for any process that signifi- 
cantly destroys the ionizing photons, unless the upper 
cut-off mass of IMF is higher than 6OM0. In Figs 6 
and 7, the tip of the vertical arrow indicates the posi- 
tion of IB model at i = if niu is increased to 12OM0. 

3.4. Classification of Sample GEHRs 

Based on the position of the GEHRs in log(-^)3 
vs {B — V)3 and {V — R)^ diagrams we classify the 
GEHRs into 5 groups. The classification is uncer- 
tain for those few regions with a poor correspondence 
between B — V and V — R colors. Nevertheless the 
classification is expected to throw more light on the 
evolutionary history of GEHRs in general. A short 
description of the five groups are given below. 

1. Young regions (tyoung < 3 Myr) with instantaneous 
burst of star formation. 

2. Moderately evolved regions (tyoung ~3-5 Myr) 
with instantaneous burst of star formation. 

3. Young regions {t < 3 Myr) with a co-existing pop- 
ulation of red supergiants from a previous generation. 

4. Moderately evolved regions (tyoung ~3-5 Myr) 
with a co-existing population of red supergiants from 
a previous generation. 

5. Evolved regions {t > 5 Myr) rich in red supergiants 
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containing stars from one or two bursts. 

Table 3 lists the objects classified under each group. 
None of the regions in our sample can be classified into 
group 5. The last row in the table gives the percentage 
of observed regions in each group. A colon in front of 
an entry suggests that the classification is uncertain 
either because of lack of correspondence between B — 
V and V — R colors or due to more than one model 
occupying a given region in our diagnostic diagrams 
(within the observational errors). The classification 
helps us to draw the following conclusions. 

1. Seventy percent of the sample regions contain stars 
from more than one burst over the last 10 Myr. The 
most recent burst has occurred within the last 3 Myr 
in 90% of these regions. 

2. The frequency of detection of young single burst 
regions is low, implying that the first generation of 
stars might be occurring deep inside the clouds, and 
hence are optically obscured when they are young and 
luminous. 

A detailed scenario of star formation, incorporating 
all the observed properties is proposed in section 5. 

3.5. Masses of Embedded Clusters 

For an assmncd IMF, the mass of the star- forming 
complexes can be estimated based on the observed 
Lyman continuum luminosity of the regions. Correc- 
tions for the decrease in ionizing luminosity due to 
evolution, destruction by dust and possible escape of 
photons have to be made in order to obtain the to- 
tal ionizing luminosity emitted from the cluster. The 
total mass (My) is thus given by, 

log(MT) = l0g((^)obs - log / - l0g{(t>)mod{t), (16) 

where {4>)mod{t) is the ionizing luminosity per unit 
mass of the cluster at an age t, as derived from the 
model for an assumed IMF. The factor / takes care of 
the loss of ionizing photons with an estimated value 
between 0.65 and 1.0. A mean ((/')mod(i) is assigned 
for each group (Tab. 3), which is used in deriving 
the mass of individual regions. We further assumed 
that only the younger burst contributes ionizing lu- 
minosity, when two populations co-exist. The masses 
derived with these assumptions and an IMF with 
a = 2.5 and m„ = 60 are given in the last column 
of Table 2. Their histogram distribution is given in 
Fig. 8 with a mean value of 

log(MT) = 5.1 ± 0.5 



for the entire sample (n = 29). While estimating 
the masses we assumed / = 1, implying complete 
detection of ionizing photons. Depending on the ex- 
act value of / for a given region, the actual masses 
might be higher by a factor of up to 2. The above 
estimates indicate that the mean mass of a cluster is 
~ 10^ Mq. The masses range by a factor of ten on 
either side of this mean value. The high mass end of 
galactic Giant Molecular Clouds (GMCs) reaches val- 
\ies up to IO^Mq (Solomon et al. 1987). From high 
resolution CO observations of molecular clouds in the 
arm of M51, Vogel et al. (1988) derive masses up to 
3 X 1O^M0, referred as Giant Molecular Associations 
(GMAs). Elmegreen (1994) has recently addressed 
the formation mechanism of such clouds and found 
that l(fMQ clouds can be formed due to gravitational 
instabilities in a spiral arm before the gas flows out of 
the spiral arm. Thus, depending on the mass of the 
parent molecular cloud, about 1-10% of the available 
mass is converted into stars in each episode of star 
formation. 

3.6. Comments on Selected GEHRs 

NGC 2363 is left out from most of the above discus- 
sions because it is different from the rest of the regions 
in the sample in terms of its low metallicity and ex- 
tinction. This region was of special interest in many 
previous investigations (Peimbert et al. 1986, Kenni- 
cutt et al. 1980). This is one of the regions with the 
smallest Balmer decrement values. The absence of 
2200 A absorption feature in the lUE spectrum (Rosa 
et al. 1985) also suggests a low interstellar extinction 
towards this region. for this region is the highest 
in the sample. The estimated age is less than 2 Myr, 
assuming an IMF with slope of 2.5 and m„ = 60. 
Broad band colors for this region are highly contami- 
nated by the strong emission lines (> 40% in V and R 
bands), and hence an investigation regarding the exis- 
tence of the second population cannot be carried out. 
Drissen et al. (1993) find the kinematic age of the ex- 
panding bubble in this region to be 2 Myr. They also 
found Wolf-Rayet (WR) stars suggesting ages around 
3 Myr. It is likely that m„ was close to 120 in this re- 
gion, with the most massive stars in their WR phase. 
The inferred mass of the cluster is 4 x IO^Mq. 

The starburst nucleus of NGC 5253 has also been 
well studied using multi- wavelength observations (Walsh 
& Roy 1989a; Gonzalez-Riestra et al. 1987; Moor- 
wood & Glass 1982). The region has sub-solar metal- 
licity. Emission lines contribute around 15-25% in 
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the BVR bands. The region agrees very weU with 
the 2-burst model, with the second population just 
formed. The inferred mass of the embedded stars is 

Among the GEHRs in spiral arms, the most inter- 
esting complex is the one lying at the tip of the north- 
east arm of NGC 2997(N2997NE-G1). The complex 
contains three knots (denoted by 1, 2, 32 in Paper I) 
of similar brightness in the continuum light while the 
knot 32 is totally absent in the Ha image. Knots 1 
and 2 are treated as independent regions and are in- 
cluded in our plots. We infer that the knot 1 itself 
contains stars from more than one burst, with a very 
young second burst. The bigger aperture measure- 
ments containing all the three knots support the re- 
sults on individual knots. A similar behaviour is also 
seen in complexes N2997NE-G3 and N2997NE-G4. 

3.7. Emission line and continuum peaks of 
GEHRs 

In two of the program galaxies (NGC 2997, 4303), 
we found a noticeable mismatch of continuum and line 
emitting peaks significantly larger than the alignment 
errors. These results are particularly worth the atten- 
tion, as there was a good agreement of the positions 
of knots in the BVR bands. As noted in Paper I, 
this mismatch has necessitated us to choose centers 
and radii of apertures by examining both the broad 
band and Ha images. We attempted getting a rough 
estimate of the separations between the emission line 
and continuum knots in two of the galaxies mentioned 
above. The measured separation between the peaks 
of nebular and stellar continuum knots is 42 ± 20 pc 
for 5 regions in NGC 2997 and 100 ± 37 pc for 13 re- 
gions in NGC 4303, respectively. It is obvious that the 
lower limit to the separations is due to the lack of res- 
olution (~ 18 and 34 pc at the distances of NGC 2997 
and 4303). The upper limit to the separation conies 
from the typical separation of star-forming complexes 
(chain of H II regions) in spiral arms. The mean sep- 
aration is 450 ± 100 pc and 1200 ± 475 pc for the two 
galaxies that we have considered. It is interesting to 
note that the separation between the neighbouring 
knots with striking brightness differences are around 
one tenth of the mean separation between distinct H ii 
complexes. Also, the separation between the nebular 
and continuum knots that we have obtained are simi- 
lar to the sizes of the largest Giant Molecular Clouds. 
This separation is also of the same order as the one 
between the old and the young populations in 30 Dor 



(Hyland et al. 1992). 

In order to explain the photometric properties of 
H II complexes, we required a population older than 
5 Myr superposed on a younger population for a ma- 
jority of the sample regions. We are tempted to be- 
lieve that in some cases we have actually resolved 
these two populations as Ha and continuum knots. 
If the older population has triggered the formation of 
younger population, we estimate the speed of prop- 
agation of the trigger as 40-100 pc in ^ 10 Myr, or 
4-10 kms~^. This value is reasonable for a trigger in 
terms of mass-loss driven stellar wind from the older 
population. 

Direct confirmation of these results will require 
spectroscopy of continuum bright knots for detectable 
signatures of red supergiants. Images taken at high 
spatial resolutions, such as possible with the Hubble 
Space Telescope will be invaluable in quantifying the 
separations with reduced uncertainties. 

3.8. Summary of Results 

From the discussions above, we draw the following 
conclusions regarding the extinction and star forma- 
tion in GEHRs in our sample. 

(1) The dereddened colors of the embedded clusters 
are found to correlate with the visual extinction es- 
timated from the Balmer decrement. Regions having 
the bluest colors have larger A„ values, implying an 
over-estimation of reddening towards embedded clus- 
ter stars. The fraction g of visual extinction experi- 
enced by the stars to that by the nebular gas is es- 
timated to be 0.7 ± 0.2, from a subset of the bluest 
regions in the sample. The value of g is likely to be 
even smaller if one considers the various uncertainties 
in the determination of g. We use a value oi g = 0.5, 
which is the value derived by Calzetti et al. (1995) for 
nuclear star-forming regions. 

(2) GEHRs with a single instantaneous burst of star 
formation are rare in our sample. Loss of ionizing 
photons amounting to ~ 35% is inferred in some re- 
gions, which may be explained as caused by the direct 
absorption by dust, if m„ = 60. Other mechanisms 
such as the escape of ionizing photons through lower 
density regions, or total absorption of some of the 
Balmer line photons by intervening dust clouds, have 
to be invoked if rn„ > 60. Escaping Lyman contin- 
uum photons possibly ionize the low density inter- 
stellar medium, thus explaining the origin of extra- 
Hii region ionized gas detected by Hunter & Gal- 
lagher (1992). 



10 



(3) There is a large scatter in the B — V and V — R col- 
ors for a given Adding a population rich in red 
supergiants, supposedly from an earlier burst, to a 
young population responsible for ionization, explains 
all the three observables. 70% of the sample regions 
are found to be in this phase of star formation. The 
estimated age of the older population is 10 Myr. 

(4) We find the Ha emitting knots spatially separated 
from the continuum knots in NGC2997 and 4303, 
from careful measurements of positions of knots on 
Ha and broad band images. If we attribute these 
knots to young and old populations, there is a spatial 
separation in addition to the age difference between 
the two populations. A mild trigger induced by mass- 
loss winds from massive stars of previous generation 
is a viable mechanism for producing successive gener- 
ations of starbursts in GEHRs. 

(5) The evaluation of IMF parameters namely m„ and 
a, which are greatly needed in understanding star for- 
mation processes in different environments, is ham- 
pered by the uncertainties in g, in quantifying the 
loss of ionizing photons and by the existence of an 
older population. The present observations do not 
compel us to choose an IMF slope far away from the 
Salpeter's value of 2.35. An upper cut-off of 60Mq 
can explain a majority of the regions. 

(6) The mass of gas converted to stars in each burst 
is 10^ Mq, with an order of magnitude spread in mass 
on either side. 

(7) The sample contains many metal rich regions. It is 

desirable that stellar evolutionary codes become avail- 
able at metallicities higher than 2Zq, considering the 
heavy dependence of evolution on metallicity. 

4. Discussion 

In this section we discuss the main results obtained 
in our work, namely the reduced extinction towards 
stellar continua and the co-existence of two popula- 
tions in a majority of the GEHRs, in the context of 
earlier work on GEHRs and related objects. There 
have been studies of GEHRs and Hii galaxies that 
show indications of reduced extinction towards stars, 
although interpreted differently by the authors. We 
look at those results afresh in the following para- 
graphs. 



4.1. H/3 equivalent widths 

The most widely used diagnostic parameter for star 
formation so far is the spectroscopically derived H/3 
equivalent width. Viallefond (1987), from his compi- 
lations, finds that the H/3 equivalent widths and the 
effective temperatures of the cluster (defined by the 
ratio of helium ionizing to hydrogen ionizing photons) 
c;annot be simultaneously explained. He finds the ob- 
served H/3 equivalent widths to be lower by a factor of 
3-4 compared to models with Salpeter IMF without 
evolution. This is essentially the same discrepancy 
that we have found in Figs 3 and 4, where observed 
log(-j^) values are much lower than the IB models 
with moderate evolution. Aperture size effects and 
usage of improper model atmosphere were thought 
by Viallefond to be partly responsible for the discrep- 
ancy. The method followed in the present study takes 
care of both these effects, but still retains the discrep- 
ancy. The data set of Terlevich et al. (1991) on Hii 
galaxies also shows the discrepancy as discussed be- 
low. 

A spectrophotometric catalogue of Hii galaxies 
containing 425 objects is published by Terlevich et 
al. (1991). It is known since the work of Sargent & 
Searle (1970) that Hii galaxies are ionized by mas- 
sive clusters of coeval OB stars. These galaxies are 
presently experiencing a high level of star formation 
activity closely resembling "starburst" galaxies on the 
higher scale and GEHRs on the lower scale. Observ- 
able properties are dominated by young stellar com- 
ponent with little or no evidence for any underlying 
population (Campbell & Terlevich 1984: Melnick et 
al. 1985). Among other quantities, Terlevich et al. 
publish Ha and H/3 equivalent widths in their cat- 
alogue. Like giant extragalactic Hii regions these 
galaxies are also thought to be simpler systems dom- 
inated by a young star-forming complex. Thus, the 
evolutionary population synthesis model discussed in 
Paper II is valid for these systems as well. As a class, 
Hii galaxies are low metallicity systems and hence 
we compare the median Ha and H/3 equivalent widths 
of this sample with our model at Z = O.IZq metallic- 
ity. The computed median Ha and H/3 equivalent 
widths for the sample regions are 172 A and 33 A 
respectively. These values compare well with mod- 
els which produce log (7^)2 values close to our me- 
dian value of —1.75. For the IMF with m„ = 60 and 
a = 2.5 at Z = O.IZq, the IB model corresponds 
to evolution beyond 5 Myr. Based on the photom- 
etry of Huchra (1977) and spectroscopy of Skillman 



11 



& Kennicutt (1993), wc derive log(-^)2 = -1.74, 
{B - V)2 = -0.02 and {V - R)2 = O.fs for the blue 
compact galaxy I Zw 18. These values suggest an in- 
stantaneous burst of star formation of age ^ 5 Myr, 
which implies that the most recent burst in IZwl8 
has properties similar to the average properties of 
Hii galaxies. Statistically one expects star-forming 
regions to be detected at their bright young phase 
(< 3 Myr). In the entire spectrophotometric cata- 
logue there are only three galaxies which have H/3 
equivalent widths > 350A, the expected value for a 
cluster younger than 3 Myr. 

With our new prescription for extinction correc- 
tion, Hq and H/3 equivalent widths depend on the 
extinction, contrary to normal convention. Correc- 
tion to this effect increases the observed equivalent 
widths, bringing them closer to the values expected 
for young clusters. The escape of Lyman contin- 
uum photons from the nebula also leads to a decrease 
in the observed equivalent widths. The addition of 
evolved stars from previous generations adds to the 
stellar continuum, thus resulting in the decrease of 
Balmer line equivalent widths. Thus, the observed 
low Balmer line equivalent widths in GEHRs and H ii 
galaxies can be explained by a combination of the 
above processes. 

4.2. UV Continuum Observations 

Young and massive stars emit most of their energy 
at UV wavelengths. Thus, the UV continuum serves 
as a very good tracer of massive star formation. The 
reduced extinction suffered by the stellar continuum, 
as we found from optical colors, should be affecting 
the observed spectra in the UV wavelength regions 
also. UV part of the spectra are essentially indepen- 
dent of IMF differences and hence are good indica- 
tors of extinction. Such spectra have become avail- 
able only recently through rocket and satellite obser- 
vations. Results from these UV spectra of GEHRs 
and other objects are discussed in the following para- 
graphs. 

Rosa et al. (1985) have published UV spectra of 
GEHRs based on the International Ultraviolet Ex- 
plorer (lUE) observations. It is found that often the 
2200 A absorption feature is weak or absent, although 
the optical extinction combined with the galactic ex- 
tinction law predicts much stronger bump. From the 
observed bumps in five objects, Rosa et al. (1985) in- 
fer E{B - V) to be not larger than 0.5 for GEHRs. 
It should be noted that the E{B — V) value obtained 



from this method measiires the "true" stellar extinc- 
tion, while that derived from the Balmer decrement 
measures the extinction towards the gaseous compo- 
nent. These observations again indicate lesser extinc- 
tion towards stellar continuum compared to the neb- 
ular radiation. 

M 83 was one of the targets for rocket UV imag- 
ing by Bohlin et al. (1990). The UV images of this 
galaxy were obtained at 1540 and 2360 A and aperture 
photometry was performed on 18 Hll regions to get 
monochromatic magnitudes mfuv and rUnuv respec- 
tively. They have performed V,B — V photometry 
based on the images of Talbot et al. (1979). By fitting 
a straight line to the reddened mnuv ~ V and mfuv — V 
colors, they find that the reddening curve towards 
stars in this galaxy is the same as that for our galax;y. 
The color excesses E{B — V) are derived by compar- 
ing the observed runuv — V colors with a model with 
Salpeter slope and m„ — 120, mi = 1.8. The mean 
value of color excess so derived is E{B — V) = 0.25 
with a small scatter (0.05 mag rms). Dufour et al. 
(1980) present optical spectrophotometry of six Hil 
regions in this galaxy apart from the nucleus. The 
derived E{B — V) values from the Balmer decrement 
have a mcidian value of 0.75, three times larger than 
the values estimated from the UV continuum. These 
observations are again consistent with our findings 
that the reddening towards stellar component is less 
than that for the nebular component. Dereddening 
the colors with the E{B — V) values of Bohlin et al. 
fits with the models with evolution within 2 Myr. 

Blue compact galaxies have attracted considerable 
attention in the UV, because of their expected blue 
spectrum. Like the study of GEHRs in M83, these 
galaxies also show much lesser extinction towards 
the UV stellar continuum compared to the values 
expected from the Balmer decrement (Fanelli et al. 
1988). 

Calzetti et al. (1994) have recently analysed lUE- 
UV and the optical spectra of 39 starburst and Blue 
Compact Galaxies in order to study the average prop- 
erties of dust extinction in extended regions of galax- 
ies. To date, this is the most exhaustive study on the 
topic of extinction towards the stellar continuum and 
nebular radiation in star-forming regions, and hence 
is the most relevant from the point of view of our 
study. Calzetti et al. find a correlation between the 
slope of the UV spectrum and the optical depth be- 
tween the Balmer emission lines Ha and H/3, which 
cannot be satisfactorily reproduced by any of the five 
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geometrical distributions of dust tlicy have consid- 
ered, for both Milky Way and LMC extinction laws. 
They derive an extinction law directly from the data 
in the UV and optical spectral range. The result- 
ing curve is characterized by an overall slope which 
is more gray compared to the Milky Way extinction 
law, and by the absence of the 2200A dust feature. 
With the new extinction law, the difference in op- 
tical depths between Ha and H/? stellar continua is 
only half of the difference between the optical depths 
for the corresponding emission lines. Our results on 
GEHRs, based on colors, are consistent with these 
results. 

Prom the above discussions it is evident that the 
star-forming complexes share the same extinction 
properties irrespective of spatial scale and strength of 
star formation or position of the star-forming site in a 
galaxy. As a rule the stellar continuum in optical and 
UV regions suffers lesser amount of extinction com- 
pared to the radiation from the surrounding nebula. 
It is likely that the same physical process controls the 
extinction in all star-forming regions. 

4.3. Thermal Radio Emission 

The extinction towards a nebula can be derived 
by using the observed intensity ratio of thermal ra- 
dio emission to Balmer lines (A^ (radio)). Most often 
the extinction is measured using two or more of the 
Balmcr linos (^4,, (Bal)). Several groups have com- 
pared the extinction derived from these two methods 
(Israel & Kennicutt 1980, Viallefond et al. 1982, Ca- 
plan & Dcharveng 1986). General conclusion from 
these studies is that Ay (radio) is around 1 mag greater 
than Ay (Bal) . Practically it is often difficult to match 
the apertures for radio and optical measurements, 
which might explain a part of the uncertainty in the 
above conclusions. However, it is believed that the 
observed higher value of Ay (radio) is real and may be 
related to the distribution of dust between the source 
and the observer. Caplan & Deharveng (1986) in- 
vestigated various dust configurations to explain ra- 
dio and optical extinction in a sample of H ii regions 
in the LMC. Prom these they conclude that a uni- 
form extinction cannot explain the LMC data. Some 
of the extinction is due to clumped dust, well out- 
side of the emission zone, and the rest is caused by 
dust located closer to or mixed with the ionized gas, 
where scattering effects are important. Measurement 
of extinction directly from resolved stars indicates 
that there is interstellar dust, quite irregularly dis- 



tributed. Higher extinction from radio measurements 
implies that some of the Balmer photons are totally 
absorbed by the intervening dust. Under these condi- 
tions, the Lyman continuum luminosity as measured 
from Balmer lines underestimates the actual ionizing 
radiation. The fraction / (or more precisely ?7esc de- 
fined in section 3.2) has a contribution from these 
undetectable Ha photons. The absolute value of / 
critically depends on the chosen m„ and hence the 
estimation of the missing number of ionizing photons 
cannot be made accurately even if one derives Ly- 
man continuum luminosities from thermal radio con- 
tinuum. 

4.4. The 30 Doradus Complex in LMC 

The Tarantula or 30 Dor nebula in the LMC is 

the closest example of a GEHR. Its relatively nearby 
location has helped its observation in a great detail. 
The near infrared observations revealed the existence 
of two populations of stars (red and blue) in this com- 
plex (McGregor & Hyland 1981). More recently Hy- 
land et al. (1992) have surveyed the region to deeper 
limits (K = 13 mag) and the earlier results arc con- 
firmed. The red population is explained to be cor- 
responding to the red supergiant phase of stars with 
mass > 2OM0 formed about 20 Myr ago. The older 
population is spatially separated from the younger 
one by distances of the order of tens of pc. Hyland et 
al. suggest that the younger cluster has formed due 
to the collision of the parent cloud with the mass loss 
wind of the older cluster stars. The present work has 
brought out the evidence for such supergiants in a 
majority of extragalactic Hii regions. 

5. Proposed Scheme of Star formation in GEHRs 

In this section, we propose a scheme of star for- 
mation, consistent with all the observed properties of 
GEHRs. The main results which the model tries to 
explain are the following. 

5.1. Observational facts 

1. The brightest GEHRs in galaxies, often preferred 

for detailed spectroscopic study, are very rarely found 
to be genuinely young with only one generation of 
stars. 

2. Among the fainter GEHRs in galaxies, there are 
a few with very high Ha emission equivalent widths 
(see Paper I). 
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3. Ha bright regions in galaxies often contain red su- 
pergiants, possibly from an earlier burst (~10 Myr 
ago) of star formation. 

4. The optical and ultraviolet continua from stars 
experience lesser amount of extinction compared to 
radiation from the ionized gas. 

5. Some of the neighbouring emitting knots have 
marked variation in the ratio of Ha to optical con- 
tinuum. 

6. On an average GEHRs contain about lO^M© of 
stellar mass in each burst. 

5.2. Proposed model 

A scenario of propagating star formation in galac- 
tic molecular clouds was discussed by Lada (1985). 
We extend this model to a larger scale to explain the 
star formation process in GEHRs. To start with, we 
summarize the known physical properties of GEHRs 
and the parent molecular clouds. GEHRs show a 
core-halo morphology in emission lines, with densi- 
ties decreasing from around 100 cm~^ in the core to 
1 cm"'^ in the halo. Kennicutt (1984) points out that 
the density in the halo almost equals the density of 
the general ISM. The average sizes of core and halo 
are around 75 pc and 250 pc, respectively (Sandage 
& Tammann 1974). In comparison. Galactic GMCs 
from which stars form have average dimensions of 
60 pc (Solomon et al. 1987), which is of the order 
of the dimension of cores of GEHRs. Interestingly, 
the Giant Molecular Associations found in the spi- 
ral arms of M51 (Rand 1993) have also dimensions 
matching that of halos, though their average densities 
are higher than the observed densities of GEHR halos. 
It is likely that GMAs contain a few GMC-like con- 
densations, surrounded by an extended low density 
inter- clump medium, which get transformed into core 
and halo respectively after massive stars are formed. 
An examination of individual regions indicate higher 
extinction towards higher Ha surface brightness re- 
gions (for example see NGC5461 in McCall et al. 
1985). Thus,the core in an Hii complex contains 
more dust than the surrounding low surface bright- 
ness halo. Emission line equivalent widths of Balmer 
lines seem to be higher near the core compared to the 
value obtained by integrating over the entire complex, 
implying that the stars are spread over a larger spatial 
extent compared to the line emitting cores of GEHRs 
(Roy et al. 1989). Stars, when born, are confined to 
the highest density regions, thus coinciding with the 



position of the Ha core. However, the relative ve- 
locity of stars with respect to the cloud make them 
diffuse out of the core into the lower density regions 
in a few million years (Leisawitz & Hauser 1988). 

With the above description of the physical proper- 
ties of GEHRs and molecular clouds, the overall pro- 
cess of star formation can be divided into four phases, 
each with a distinct observable signature. These four 
phases are schematically depicted in Fig. 9. The outer 
boundary denotes the parent molecular cloud, whose 
dimension is around 300 pc. The GEHR core is shown 
by the circle, with the bigger concentric circle denot- 
ing the extent of the ionized halo. Stars are denoted 
by asterix. The irregular pattern indicates the den- 
sity inhomogeneities in the parent cloud. The approx- 
imate time period in each phase is indicated. 

Phase 1: We assume the GEHRs to form in Giant 
Molecular Associations, which contain several regions 
of higher than average densities. Some of these con- 
densations grow into GMC-like structures, the dens- 
est of which collapses to form stars. This results in 
an ionization of the high density gas in the imme- 
diate vicinity, forming the core of the GEHR. Dust 
associated with the molecular cloud offers consider- 
able optical extinction, allowing only a few of the 
young regions to be detected optically as faint Ha 
emitting knots with weak or no continuum. Regions 
with high Ha equivalent widths are in this phase of 
evolution. The large extinction when the cluster is 
young, makes its optical detection very difficult, 
is expected to be equal to Tg in the very early phase, 
with a tendency for to decrease as some of the op- 
tically luminous stars diffuse to low density regions. 
The ionization front is advancing into low density re- 
gions. This phase might last for 3 Myr at the most. 
Phase 2: There might be a net reduction in the 
extinction towards a GEHR in a few million years, 
as the energetic events associated with star forma- 
tion destroy the dust particles in the core. The re- 
duced extinction makes these regions optically visible 
for the first time. Embedded cluster stars might expe- 
rience additional reduction in extinction as more and 
more stars diffuse out of the high density cores. By 
this time the ionization front has also advanced past 
the high density regions of the molecular cloud, thus 
creating an extended low density ionized halo, which 
can be detected only through sensitive observations. 
This might lead to an underestimation of the ioniz- 
ing photons. However, emission from the core still 
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dominates the total Balrncr emission from the region. 
Under this configuration, photons from ionized gas 
encounter more dust clumps, compared to stellar pho- 
tons. Thus, the stellar continuum experiences lesser 
extinction compared to radiation from the ionized gas 
(rg < Tg). This phase might continue up to around 
5 Myr. The second sketch in Fig. 9 depicts this stage. 
Phase 3: The onset of this phase is characterized 
by the appearance of luminous red and/or blue super- 
giant stars in the cluster, as the massive stars evolve 
off the main sequence in around 5 Myr. The cluster 
appears brightest in the stellar continua from ultravi- 
olet to near-infrared during this phase. This luminous 
period continues until all the massive stars go through 
the red supergiant phase (~12 Myr). In this phase the 
region is very faint in Ha. 

Massive stars lose a large amormt of mass and en- 
ergy during their life time, both on and off the main 
sequence. This results in an increased pressure on 
the gas left behind after the initial burst of star for- 
mation, thus creating another density condensation 
in the same GMA. Evidence is accumulating for the 
presence of left over molecular gas in the vicinity of 
GEHRs {e.g. NGC604 in M33: Israel et al. 1990; 
Wilson & Scoville 1992). The shocked layer starts col- 
lapsing, eventually forming a new generation of stars. 
The delay between the onset of collapse and genera- 
tion of new stars is given by (Elmegreen 1992), 

.^.0.25 

fdelay = 2.4 — j Myr. (17) 

For a cloud of mass lO^M0participating in the col- 
lapse this corresponds to ideiay = 4.3 Myr. The arrow 
in the third panel indicates that the cloud is collaps- 
ing. 

Phase 4: The end of phase 3 and the beginning of 
phase 4 are characterized by the birth of the second 
generation of massive stars the same parent molecular 
cloud. Typically this happens at around 8 Myr, with 
the exact epoch depending on the mass of the parent 
cloud and the onset of collapse. During this phase, 
the region becomes brighter again in Ha. By this 
time the stars from the first generation have evolved, 
and the cluster itself has moved into regions of lesser 
extinction. In this phase, the bulk of observed stel- 
lar continuum originates from the older generation 
and hence suffers lesser extinction compared to the 
ionized gas situated in dense regions, thus continuing 
the trend shown in Phase 2. GEHRs are expected to 
be in this phase at about 10 Myr since the beginning 



of the first generation of stars. For typical propaga- 
tion speeds (~ 10 kms~^), the separation between the 
two populations amounts to ~ 100 pc. Thus the two 
generations of stars can be distinguished, under good 
spatial resolutions, as neighbouring Ha and contin- 
uum (stellar) peaks. 

In this study, we restricted the evolution of GEHRs 
up to the second generation. In reality, star forma- 
tion might continue until the parent clouds of GEHRs 
are fully disrupted by the energetic events such as 
multiple supernova explosions. It is likely that big- 
ger systems like Hll galaxies and central regions of 
galaxies experience several bursts of star formation 
before the consumption of gas is complete. However, 
the process of formation of multiple bursts might be 
different in these bigger systems. Recently, Korcha- 
gin et al. (1995) have constructed propagating star 
formation models for the nuclear hot-spots in galax- 
ies and found the star formation lasting for a few tens 
of million years. 

Before building on the above tentative model, it is 
important to establish the different phases by elab- 
orate observations. This would involve a detailed 
study of selected high Balmer equivalent width re- 
gions to establish that they are young and also that 
their stellar continuum experiences similar extinction 
as the radiation from the ionized gas. Regions in the 
present sample with definite indication of older gen- 
eration should be observed for the signatures of red 
supergiants. GEHRs should be imaged in the opti- 
cal emission and continuum bands with high spatial 
resolutions, in order to directly resolve the spatially 
separated generations of stars, one bright in the emis- 
sion line and the other bright in the stellar continuum. 
As all these observations are well within the observ- 
ing capabilities of most medium-sized optical-infrared 
telescopes, we hope that the above model will be put 
to test soon. 

6. Conclusions 

We find an evidence for reduced extinction to- 
wards stellar continuum compared to the extinction 
derived from the Balmer decrement. Genuinely young 
GEHRs containing stars from only one generation 
seem to be very rare in the present sample. Instead, 
the young and moderately evolved regions contain an 
older population rich in red supergiants. All these ob- 
servational facts are put together resulting in a simple 
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qualitative model, which can be tested by a detailed 
study of selected regions. In the model, GEHRs are 
obscured by the dust associated with the high den- 
sity gas when they are young, and arc detected only 
when the stars have moved out of the highest density 
regions. This leads to a lower effective optical depth 
towards the stellar continua compared to the emis- 
sion lines from the ionized gas, whenever GEHRs are 
detected. The mass-loss wind from the massive stars 
triggers a second burst of star formation in the parent 
cloud, providing a fresh input of ionizing radiation. 
The presence of stars from two bursts and reduced 
net extinction makes these regions easily noticeable 
on galaxy images, and hence are preferred regions in 
most of the studies. Star formation in GEHRs stops 
after two or more bursts, probably because of disrup- 
tion of the parent gas cloud. 

The suggestions made by the referee Marshall L. 
McCall has helped enormously in improving the pa- 
per, especially in revising the proposed model for star- 
formation. 
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Fig. 1. - Observed GEHRs in log(-^)i vs {B - 
V)i plane. The nebular contribution is subtracted 
from broad bands and extinction corrections are done 
only for galactic dust. Typical rms errors are shown 
by the cross at the bottom left corner. The arrow 
indicates extinction correction corresponding to Ay = 
1 mag. The evolution of clusters at Z = 0.1, 1 and 
2Zq are shown by the short-dashed, long-dashed and 
solid curves respectively. An IMF with parameters 
niu — 60 and a = 2.5 is used. The dots on the line 
are spaced 0.5 Myr apart with the dot corresponding 
to 4 Myr marked by the letter 4. 

Fig. 2. — Observed regions in log(;^)i vs {V — R)i 
plane. See Fig. 1 for other details. 

Fig. 3. — Observed regions in log(-^)2 vs {B — V)2 
plane. Extinction corrections are done based on the 
Balmer decrement. Other details are as in Fig. 1, 
except that only Z = 2Zq model is shown and the 
evolution is shown only upto 6 Myr. 

Fig. 4. — Observed regions in log(-^)2 vs {V — R)2 
plane. See Fig. 3 for other details. 

Fig. 5. — The B — y colors after correction for extinc- 
tion derived from the Balmer decrement are plotted 
against Ay corrected for galactic extinction. Note the 
trend of bluer regions having larger Ay. 

Fig. 6. — Observed regions in log(-^)3 vs {B — V)3 
plane. Extinction corrections are done using Aye 
with g = 0.50. Solid line represents the Z = 2Zq 
model evolved upto 6 Myr. The same model is shifted 
by 0.20 dex (see text for details) and shown by the 
dashed line. A model with rn„ = 120 and no shift 
begins at the tip of the arrow. Sequences with short- 
dashed lines correspond to composite models with the 
younger population at (top), 3 (middle) and 4.5 
(bottom) Myr age superposed on an older popula- 
tion. The older population is 6-14 Myr old along the 
sequence. Each sequence begins with a dot marked 
by vertical line. 

Fig. 7. — Observed regions in log(-^)3 vs {V — R)^ 
plane. See Fig. 6 for other details. 

Fig. 8. — Histogram of masses of GEHRs. 

Fig. 9. — A schematic diagram of the four phases 
during star-formation in GEHRs. The outer bound- 
ary in each sketch represents a Giant Molecular As- 
sociation which contains a GEHR, surrounded by the 



low density gas. The GEHR itself contains cluster 
stars, clumpy gas and dust as shown. The formation 
of GEHR core, subsequent formation of the halo, the 
spreading out of optically luminous stars into lower 
density regions, and subsequent triggering of a second 
burst of star formation in the left out gas are some of 
the salient features of the model. Typical time scales 
involved are marked on the top right. See section 5.2 
for details. 
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